Two interacting electrons in a magnetic field: 
comparison of semiclassical, quantum, and 
variational solutions 

Tobias Kramer 

Institutfiir Theoretische Physik, Universitat Regensburg, 93040 Regensburg, Germany 

Abstract. The quantum mechanical many-body problem is rarely analytically solvable. One no- 
table exception is the case of two electrons interacting via a Coulomb potential in a uniform mag- 
netic field. The motion is confined to a two-dimensional plane, which is commonly the case in 
nanodevices. We compare the exact solution with the semiclassical energy spectrum and study the 
time-dependent dynamics of the system using the time-dependent variational principle. 
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INTRODUCTION 

The study of the motion of charged particles in an external magnetic field is of impor- 
tance in many fields of physics. Astrophysical simulations need to model plasmas, while 
in solid-state physics at low temperatures magnetic fields cause the integer and fractional 
quantum Hall effects. In all cases the understanding of interactions is required to de- 
scribe electronic transport since the motion of charged particles in an external magnetic 
field is driven by the electric fields generated by the charges themselves, and applied 
externally [1, 2]. Few-electron systems can be engineered with great precision in nan- 
odevices in quantum-dots and the excitation spectrum of the system has been studied 
experimentally and theoretically [3]. The presence of interactions causes intricate cor- 
relation effects, most prominently discussed within the fractional quantum Hall effect. 
But before one can proceed to large systems, it is important to reduce the complexity 
and study the simplest example, which still permits an analytical and numerical exact 
solution. For two- and many electron-atoms Moshinsky, Novaro, and Calles devised 
the pseudo-atom [4, 5], where the Coulomb interactions are replaced by a repelling har- 
monic oscillator potential in order to solve the many-body problem exactly and to obtain 
quantitative error bounds for the Hartree-Fock method. Interacting electrons in a mag- 
netic field show a particular rich behavior, since the repulsion leads to an initial increase 
in distance. This is converted by the velocity-dependent Lorentz force into the cyclotron 
motion, which brings all charges back together and no particle escape is possible if the 
system is confined to two-dimensions. The solution of Fock [6], Landau [7], and Darwin 
[8] for the limit of the non-interacting case are discussed in many textbooks on nanode- 
vices [9, 3], while the solution for the interacting case has been found many decades 
later [10, 11]. Already the classical system reveals a complicated dynamics due to the 
interplay of the Lorentz force in the magnetic field and the electric field generated by the 



electron-electron repulsion [12]. In nanodevices, electronic properties are often probed 
in a two-dimensional layer perpendicular to the magnetic field direction [9, 13] and thus 
a reduction of the three-dimensional problem to a quasi two-dimensional one is justified 
and needed to understand the experimental measurements. 

Here, we analyze the interacting two-electron problem in the absence of an addi- 
tional confining potential, since we are interested in the intrinsic dynamical properties 
of the system. After studying the exact quantum-mechanical stationary solution, which 
is available only for specific parameters, we present a semiclassical analysis covering 
all interaction values. Then we present a time-dependent, dynamical view of the elec- 
tron dynamics in the frame-work of the time-dependent variational principle (TDVP) 
and compare the TDVP to the numerical solution. The two-body part of an interacting 
quantum-mechanical systems already displays profound differences to the classical dy- 
namics caused by the antisymmetrization requirement for fermionic dynamics [14]. Fi- 
nally we present an outlook how to proceed from the present two-body analysis to more 
particles in a systematic fashion, which is required to study transport in nanodevices. 



THE SPECTRUM OF TWO INTERACTING ELECTRONS 

Ignoring the extension of the wave-function along the third dimension in the direction 
of the magnetic field, the Hamiltonian for the two electrons is given by 

with Larmor frequency Ol — ^ and angular momentum L^^i — —ih{xidy^i — yidxj), 
i = 1,2. The strength of the Coulomb interaction is given by the parameter Jc. The 
Hamiltonian Hee separates into a center-of-mass term Hem and a relative term flrel in 
Jacobi coordinates 

p2 J 

^cm = — + -M(ofR^-(OiLz, 

^rel = 24l^2^^^^^~^^'^]^y 

with respective coordinates and masses 

r = r2-ri, (4) 

P = -i^Vr, (5) 
ju = m/2. (6) 

Only the relative part of the Hamiltonian contains the Coulomb interaction, and the 

solutions for the non-interacting center-of-mass Hamiltonian Hem are Landau levels with 
eigenfunctions and energies given by products of exponential functions and Hermite 
polynomials [15]. The ground state attains a gaussian shape with the eigenfunction and 
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eigenenergy given by 

W2(X,i') = ^exp(-^^), A = ^, E^=m. (7, 



Exact solutions for particular interaction strengths 



The solutions of the relative part H^ei are available in closed form for specific ratios 
of the Coulomb interaction and the magnetic field strength [10, 11, 16]. Only in these 

cases the eigenstates are given by the product of a polynomial of finite degree and a 
decaying exponential function. However we do not view the absence of finite-degree 
polynomial solutions as an indication for the appearance of anyons (as done in [10]), 
but find numerically and semiclassicaly normalizable solutions for all values of the 
interaction strength. The conditions for a closed form solution are related to a hidden 
symmetry within the problem [17, 18]. Here, we list the first solutions classified by the 
angular momentum / and the quantum number rir 
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(jc,};) = A^^"-') i?(«'-)(r)e^^',r= V^2^/, 8 = arctan(A;,3;), / = 0,±1,±2, . . . (8) 



up to the normalization factor Ar("'^'0 . We set ^ = /i = ft)/ = 1 in Eq. 3 and thus all energies 
are given in units of h(Oi, all lengths in units of 2W ^, and the Coulomb interaction 



strength in units of ^ \/ x 
For K = v'|/| + l/2 
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for / < 0, 
for / > 0. 



Semiclassical WKB spectrum 



To construct the entire energy spectrum at interaction values between the ones tab- 
ulated above requires numerical methods or alternatively to employ approximation 
schemes. The semiclassical energy spectrum can be obtained for all interaction strengths 
by following Ref. [19], where we again set the external confinement to zero. The ansatz 
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leads to the radial equation 

X"ir) + 
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(13) 



Upon making the Langer correction appropriate for two-dimensions, — 1/4 — > [20], 
the semiclassical quantization condition reads 



dr\ 2E + 21 



— -r2 = ;?r(n+l/2), 
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(14) 



where ri^2 denote the inner and outer turning points (zeroes) of the radial equation. The 
polynomial in the square root has four zeroes, a, b, c and d, which are ordered such that 
b> a and c, J < 0. For given k, n and m we obtain an implicit equation for the energies 
Once a, b, c, d are known, the integral (14) is given in Ref. [19] 
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and written in terms of the complete elliptic integrals K, IT 
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The semiclassical spectrum as a function of k is shown in Fig 1. The ellipsoids mark 
exact results from the quantum solution at specific values of Jc. The dashed lines denote 
a reversed sense of rotation (/ < 0) compared to the direction preferred by the Lorentz 
force. The states with a reversed sense of rotation / < are shifted to higher energies. Our 
results are in line with the Landau level degeneracy for the non-interacting case, where 
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FIGURE 1. (Color online) Semiclassical spectrum E{k). Color denotes the value of the angular mo- 
mentum / (solid lines I > 0, dashed lines / < 0, black: I — 0, blue: / = ±1, green I = ±2, red / = ±3, 
magenta I = ±4, yellow I — ±5. The ellipsoids denote exact quantum-mechanical solutions based on 
Heun polynomials. 



one obtains [15], p. 193 that the Landau level of quantum number n (r = 0, 1 , 2, . . . , oo) 
contains the angular momenta Z = — «, — n + 1 , —n + 2, . . . , oo. The Coulomb interaction 
lifts this degeneracy and reveals the number of negative angular momenta states con- 
tained in each Landau level. Note that the Coulomb interaction K depends in Sl-units on 
the magnetic field. The spectrum for fixed Coulomb interaction as a function of magnetic 
field is show in Fig. 2. 

TIME-DEPENDENT SOLUTION 

In order to understand the role of interactions and antisymmetrization in scattering 
effects between two electrons, it is instructive to consider the time-dependent dynamics 
of two initially separated electrons. For simplicity we restrict ourselves to cases where 
the center-of-mass part of the system remains in the ground state and we only consider 
the time-evolution of the relative part. As initial orbitals for the two electrons, we pick 
two shifted Gaussian functions which provide an exact solution in the absence of the 



5 10 15 20 

magnetic field [Tesla] 

FIGURE 2. (Color online) Semiclassical spectrum E (B) / (hcoi) as function of magnetic field B in Tesla, 
(effective electron mass m* =0.1, dielectric constant e = 10). Color denotes the value of the angular 
momentum I (solid lines / > 0, dashed lines / < 0, black: / = 0, blue: / = ±1, green / = ±2, red / = ±3, 
magenta / = ±4, yellow / = ±5. 




Coulomb interaction 

^ -l^expl-^l^^l. (21, 

From the two orbitals we form the symmetric and antisymmetric and not- symmetrized 
combinations for the two-electron wave function 

i/(ri,r2) = [r(ri)(^)^(r2)-r(r2)(^)^(ri)]/v^ (22) 
V/^(ri,r2) = [r(ri)(^)^(r2)-fr(r2)(^)^(ri)]/v^ (23) 
V/^(ri,r2) = [r{ri)^Hr2)]. (24) 

Since the orbitals are not orthogonal to each other, the normalization N differs from 
\/2. The width parameter a of the Gaussian is chosen such that the electrons are in the 

lowest Landau level a = ^J^- We place the two electrons initially at rest (this requires 
a gauge factor!) at a separation ofd = 2xo- Next we transform the states to relative and 



center-of-mass coordinates, where the wavefunction becomes a product 

V^(ri,r2) = V^rei(r2-ri)v^cm((ri+r2)/2). (25) 
In relatives coordinates we obtain (ii = h = coi — I) 



V4i(r) = J-^—smh{xxo)c-^^+y'y\ (26) 
K{e^o — 1) 



V'^i(r) = J —-r—-cosh{xxo)c-(^'+y'^/\ (27) 
y ;r(e-^o + 1) 

V^i(r) = y|e-«--")'+^')/2, (28) 
and the center-of-mass wavefunction reads 

WS(X,F) = Ae-2(^^+^^). (29) 

The total wave-function is then recovered via 

W^^''^^{xi,yi;x2,y2) = ¥^i{xi+X2)/2,{yi+y2)/2)x(rff'^{x2-xi,y2-yi). (30) 
As time-dependent observable quantities we track the density of the relative motion 

Prei(r;0= |v^rei(r;OF (31) 

and the single-particle density obtained by integrating over one set of particle coordi- 
nates 

Pspd(ri;0 = ydr2|V'(ri,r2)|^ (32) 
= y"dr2|v^rei(r2-ri;0lVcm((ri+r2)/2;0p (33) 
The single-particle densities of the initial states are 

Ps1,d(r) = ^^^2_^^ (e^''/"cosh(2x^)-l) (34) 

Ps1,d(r) = , o ^. (e^/"cosh(2x^) + l) (35) 
^ ;r(e^''-|-l) ^ / 

2e-2{{x-Xo/2f+y^) 

Pspd(r) = , (36) 

with normalization 

/drPst,f'^(r) = l. (37) 



Numerical time-propagation 



The initial state in the relative coordinate can be represented in discretized form on a 
regular two-dimensional mesh, and the time evolution is performed using a split-operator 
method adopted to the presence of a magnetic field [21]. The Hamiltonian is represented 
in a mixed momentum-position representation necessary due to the angular-momentum 
operator L^. For small time-steps A? the time-evolution is given by 

r '\N 
\l/{r,t' + NAt) fa Q-iAt/nV{xj)/2 ^-l^~iAt/hTpy^^^j^-l^-iAt/nTpj,^y j^^^~iAt/nV{x,y)\ 

xe'^/^^(^'^)/V(r,^0, (38) 
where denote partial Fourier transforms with respect to only one-dimension and 

2 2 

7>.,y=^-®LPxy, 7>„x=0 + «LPyX, V(x,3;) = ^M®/V+^. (39) 

The coherent- state propagator for vanishing Coulomb interaction is known in closed 
form [22] and can be used to test the numerical scheme. 



Time-dependent variational principle 

Next we consider the dynamics of two electrons based on the application of the 
time-dependent variational principle (TDVP) due to Dirac and Frenkel, which has been 
reviewed in [23, 24]. The starting point is the action 
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The principle of extremal action requires dS — J dt dL = 0. For the Gaussian orbitals 
used in the previous section is is convenient to use a notation related to coherent states 
and we label the^arameters of the variational wavefunction by z: \i/= ¥{z) = \z)- Using 

idt = 5 (- + ) we get 



5L = 







z 


z)-{z 


H 


z) 


{z 


z) 



■{z\z) - -{z\z) - {z\H\z) 



{8z\z) 
{z\z)^ 



(41) 



-l-c.c. 



The time-derivatives of the form 5i are removed using the identity 
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and we rewrite 5L 
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{5z\z) {5z\H\z) {5z\z) 
{z\z) {z\z) {z\zf 



[i{z\z)-{z\H\z)] +C.C., 



(43) 



where we dropped all total time derivatives since the variation vanishes at ti and t2. 
Finally we obtain from the variation with respect to dz 
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(44) 



the equations of motion for the parameterized wave function by denoting 

{W{z)\H\xif{z)) 



|V/) = |V/(zi...zyv)), N^{w{z)\w{z)), ^ = 



{V{z)\v{z)) 



(45) 



and introducing dz = T^id^jSz, 
5L = £ 



5zl (46) 



For a compact notation we introduce the matrix 

<92 



InN, 



(47) 



which allows us to recast the variational equations in a shape familiar from Hamilton's 
equations 
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In the ansatz for the variational function the phase and normalization were explicitly 
removed. A normalized state with phase is restored by multiplying the variational state 
xi/hy an exponential factor 
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Choice of variational wave-function 

Any variational principle can return only good results if the variational functions are 
chosen with physical insight into the problem at hand. The choice of coherent states for 
electron-orbitals in a magnetic field is physically motivated by their property to restore 
the exact result for the non-interacting case. The Coulomb interaction distorts the coher- 
ent state evolution by changing the trajectories and momenta encoded in qx,1y, Pxi Py 
In addition the width of the wave-function is no-longer fixed and we generalize our 



"frozen" coherent-state variational functions to include squeezed coherent states with 
adaptable "thawed" widths. A variable width requires to consider two more variables, 
since the time-evolution leads to a dispersion which has to be represented by a complex- 
valued parameter [14]. In total we describe the relative motion by six variational param- 
eters qxi qy,Px, PyiU,w based on the function 

wg,{x^y;p..Py^q.Ay.u.w) = e-^((^-^)'+(^-^)')+^^^^+^^^^ (50) 
which we antisymmetrize 

Wtiix,y;px,py,qx,qy,u,w) = \i/^^i(x,y;px,Py,qx,qy,u,w) 

-v^iC^^y; -Px, -py, -qx, -qy, -w, -w) (5i) 

in order to describe two electrons with the spin pointing in the same direction. It is 
instructive to compare the expectation value of the Hamiltonian with respect to i/^j and 
V^^j. In the first case we obtain 

pI + pI 1 9 9 A/^e~*^/('^'^^) 9 9 
+ r - (OiiPyqx-Pxqy) + ^ 7^ Io(r/(4w^)) 
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"classical" 



^ 1 2.,^J2 



+ — + —^+w^^cof, (52) 

" V ' 

"uncertainty" 

where the the first parts corresponds to the Hamiltonian of the classical charge den- 

sity |i//rei| and the last part adds the minimum uncertainty and represents a quantum- 
mechanical effect. The antisymmetrized expectation value contains in addition correla- 
tion and exchange terms 

= '<Vl"'°'y' (53) 

_ Px+Py 1 ,,-2 2 / \ I ^ I 2, ,,.2 



+k: 1 + 



+ ^ MO/ r - 0)1 {Pyqx - Pxqy) + t^tjTT + 

1 



„2 //•/i,.,2 



^-1 



e 



-Io(?7(4w^)) 



where ^ denotes a distance in phase space 
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^ = 2 [{uqx - wpxf + {uqy - wpyf^ + ^""^2^ ■ (54) 



For large values of ^ the antisymmetrized expectation value of the Hamiltonian ap- 
proaches J^'^. In order to complete the derivation of the equations of motion, the matrix 
C has to be determined. The resulting expressions are lengthy and not reproduced here. 
For the antisymmetrized case, all variational parameters are now coupled, whereas for 
the unsymmetrized matrix a fully sympletic relation holds for the position and momenta 
as well as for the width parameters 
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As for the expectation value of the Hamiltonian, the antisymmetrized C matrix ap- 
proaches the symplectic form for large distances in phase space. The system of cou- 
pled differential equations for the time-evolution of the initial state is now completely 
determined 
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but the integration has to be carried out numerically. We compare two different cases to 
the exact propagation in order to determine the quality of the variational function. The 
first case restricts the variational degrees of freedom to only four parameters, by keeping 
w fixed and setting m = 0, while the second approach includes all six parameters. The two 
TDVP solutions are compared to the numerically exact solution which we obtain with 
the wave-packet propagation scheme described above. In Fig. 3 we show the density of 
the relative wave-function for different propagation times. After starting from the same 
initial state, the TDVP solutions deviate with longer times. The frequency of rotation 
around the origin is a sensitive measure of the shape of the wave-function, since a slight 
change in the width of the wavefunction alters the electric field and affects the Lorentz- 
force induced rotation. A mechanical analog is the change of rotation rate of a spinning 
object due to changes of its moment of inertia. The thawed (squeezed) coherent-state can 
adapt its rotation rate better than the frozen- width coherent-state and thus stays closer to 
the exact solution for a longer time-evolution. 



OUTLOOK: FROM TWO TO MANY ELECTRONS 

The two-electrons is a special case of the general many-body problem. The addition of 

more particles complicates the problem considerably and no exact analytic solutions are 
known for the general case. In the case of a classical system of interacting particles, the 
self-generated fields of the electronic charges are the foundation of the classical Hall 
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FIGURE 3. Comparison of the time-evolved densities li^reK-^jy)! ^^^J" ^ different number of prop- 
agation steps (time step At = 0.02, Coulomb interaction k = I, initial offset Xo = l/\/2, initial width 
w = 1/ \/2). The left column shows the exact result, the middle column the "thawed" wavepacket of a 
squeezed coherent state with all six variational degrees of freedom, and the right column the fixed-width 
("frozen") TDVP result for a coherent state. With increasing time the TDVP solutions start to deviate from 
the exact solution. 



effect. We have demonstrated in [1] that realistic boundary conditions and the Coulomb 
interactions completely alter the current-density distribution and produce the Hall field, 
in contrast to the non-interacting solution. The inclusion of quantum-effects in a classical 
calculation requires to account for the Pauli-principle and to recover the Fermi-Dirac 
statistics of fermionic systems. The TDVP offers a possibility to systematically move in 
this direction [14] which we are currently exploring. 
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